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Abstract 

We analyse the centrality dependence of thermal parameters deduced from hadron mul- 
tiplicities in Au + Au collisions at ^/snn = 130 GeV. While the chemical freeze-out 
temperature and chemical potentials are found to be roughly centrality-independent, the 
strangeness saturation factor 75 increases with participant number towards unity, sup- 
porting the assumption of equilibrium freeze-out conditions in central collisions. 
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I. INTRODUCTION 

Statistical-thermal models (cf. Q| for recent surveys and references therein) have 
enjoyed remarkable success in describing hadron multiplicities observed in both 
heavy- ion and elementary collisions over a wide range of energies 
The final state multiplicities are reproduced in these models with very few param- 
eters. The prominent ones are the chemical freeze-out temperature T and baryon 
chemical potential fiB- There are further chemical potentials to be considered be- 
low. A compilation of the freeze-out parameters over the accessible range of energies 
reveals a continuous curve in the T — plane which may be characterized by an 
energy per hadron of the order of 1 GeV [4|]. While the occurrence of such a uni- 
versal freeze-out curve is interesting in itself and is useful for extrapolations and 
interpolations as well, its particular meaning could be the approximate coincidence 
with the confinement border line, at least at high energies, as conjectured in 

We do not touch here on the disputed question as to why also the multiplicities 
of the rare multi-strange hadrons seem to obey chemical equilibrium (see jsf). In- 
stead, we focus on one possible indicator of incomplete equilibrium, the strangeness 
saturation factor 7^. This factor has been introduced H to account for an appar- 
ent under-saturation of strange hadrons. Various previous analyses ^Ifl] considered 
such a factor as necessary to accomplish a satisfactory description of data. For SPS 
energies, ^Jsnn = 6 ■ ■ ■ 20 GeV ^ it has been shown that 7s grows with increas- 
ing beam energy and centrality. First attempts based on very restricted data 
sets, find also at RHIC energies a preliminary indication of rising 7^ with increasing 
centrality. It is the subject of this paper to investigate in detail the dependence 
of the thermal parameters including 7^ on the centrality in collisions of Au+Au at 
y/SNN = 130 GeV. By now there is a sufficiently large basis of published data at 
this energy to allow a thorough analysis. 

Our paper is organized as follows. In section II we recall the formulation of our 
employed statistical-thermal model. The analysed data are listed in section III. Our 
fit procedures are described in section IV and results are presented in section V. 
The summary and conclusions can be found in section VI. 
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II. STATISTICAL-THERMAL MODEL 



For suitably large systems at RHIC collider energies, hadron multiplicities 
are analysed by employing the grand-canonical partition function Z{V,T, fii) = 
Tr exp{ — {H — J2i f^iQi)/T} , where H is the statistical operator of the system, 
T denotes the temperature, and /ij and Qi are respectively the chemical potentials 
and corresponding conserved charges in the system. To be specific, we use here 
the baryon charge (corresponding to fis) and the strangeness charge (corresponding 
to fis) as unconstrained quantities. (There are other possible chemical potentials, 
e.g., /ig as electric charge potential. However, we put /ig = with the reasoning 
explained below.) 

The primordial particle numbers are accordingly. 



BW(mi), (2.1) 

where we include phenomenologically the strangeness saturation factor 7^ (with \Si\ 
the total number of strange and anti-strange quarks in hadron species i) to account 



for the possibility of incomplete equilibration in this sector; Ei = yp^ + m? and 
BW(mj) is the Breit-Wigner distribution (to be replaced by a 5-function for stable 
hadrons) to be integrated from thresholds with appropriate widths Fj (in practice 
the interval [max {threshold, mj— 2Fj}, mj+2Fj] is sufficient). The fiducial volume V 
drops in particle ratios, gi is the degeneracy factor of particle species i with vacuum 
mass rrii. (The study of the effect of in-medium masses (cf. |13]) is interesting 
because it might shed light on the chiral properties of the medium created in heavy- 
ion collisions at RHIC.) Eq. (j2.H) simplifies in the Boltzmann approximation (i.e., 
discarding the spin-statistics factor ±1) and when neglecting the energy distribution 
of resonances. We use here the full expression Eq. ()2.1|) and include all hadron states 



with u,d,s quarks and corresponding anti-quarks up to 1.7 or 2.6 GeV 32] with 
masses and total widths according to the particle data group listing The final 
particle numbers to be compared with experiment are Ni = Nf"'^ + J2j Bi^^^Nj"'^ 
due to decays of unstable particles with branching ratios Br-'^'. 

Originally, such a description was thought to be justified for multiplicities mea- 
sured over the entire phase-space, since many dynamical effects cancel out in ratios 



of hadron yields 



At sufficiently high energy, such as at RHIC, however, ratios 



of mid-rapidity yields are also found to be well described by the statistical-thermal 







model m a 02 
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III. THE DATA 

In contrast to SPS energies, where both mid-rapidity and fully phase space in- 
tegrated data are at our disposal (see 3| for a comparison), at RHIC energy of 
y/SNN = 130 GeV sufficiently many data of identified hadrons are available only 
at mid-rapidity. The data analysed here were compiled from PHENIX, PHOBOS, 
BRAHMS and STAR experimental results. These experiments differ in their cuts 
on the data. It turns out that the results for the deduced thermal freeze-out pa- 
rameters depend in some cases sensitively on the employed data sets. Therefore, for 
definiteness Table E lists the ratios used as a function of collision centrality after re- 
calculation to a common centrality binning as in Subscript "(1)" denotes yields 
corrected for weak feed-down, while yields labelled by "(2)" include a contribution 
from weak decays. Three centrality bins were selected with participant numbers 
63.5 ± 8.4 (peripheral), 210.5 ± 8.4 (mid-central) and 317.0 ± 8.2 (central). 

IV. THE FIT PROCEDURES 

According to our propositions, T, /i^, /ig and 7^ are unconstrained fit parameters. 
For each centrality class, four fits were performed: 
Fit I: all available ratios listed in Table I included; 
Fit II: only ratios of 7c, K, p, A and S included; 
Fit III: only ratios of tt, K, p and A included; 
Fit IV: only ratios of tt, K and p included. 

In this way, the effect of the various multiplicities on the thermal parameters are 
investigated. (Various other fitting procedures are considered in Q|.) 

The final ratios compared with experiment include both a primordial and a 
decay contribution. As mentioned above, in some cases the experimental data have 
been corrected for weak decays. However, where such corrections have not been 
made, the infiuence of weak decays is included in the following way: 
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TABLE I: Ratios of hadrons in Au + Au collisions at y/sj^N = 130 GeV from various 
RHIC experiments. 



Ratio 


Experiment 


Central 


Mid-Central 


Peripheral 




BRAHMS [19] 


0.990±0.100 








PHENIX [20, 211 


0.960±0.177 


0.920±0.170 


0.933±0.172 




PHOBOS [22] 


1.000±0.022 








STAR [17] 


1.000±0.073 


1.000±0.073 


1.000 ± 0.073 


^{2)1-^(2) 


PHENIX [2ii, 
PHOBOS [22] 


1.152±0.240 
1.099±0.111 


1.292±0.268 


1.322±0.284 




STAR [23] 


1.109±0.022 


1.105±0.036 


1.120±0.040 




PHENIX [20, 21] 


0.680±0.149 


0.671±0.142 


0.717±0.157 


F{2)l F(2) 


BRAHMS [19] 
PHOBOS [22] 


0.650±0.092 
0.600±0.072 








STAR [26] 


0.714±0.050 


0.724±0.050 


0.764±0.053 


Wl w 


PHENIX [2^1 


0.750±0.180 


0.798±0.197 


0.795±0.197 


[2)1 [2) 


STAR [18] 


0.719±0.090 


0.739±0.092 


0.744±0.100 




STAR [28,^ 


0.840±0.053 


0.822±0.114 


0.815±0.096 




STAR [29, 30] 


1.062±0.410 






{2)1 (2) 


PHENIX [2ii, 21] 


0.151±0.030 


0.134±0.027 


0.116±0.023 




STAR \XL2^ 


0.151±0.022 


0.147±0.022 


0.130±0.019 


K%l tt7^, 

bl (2) 


STAR [17:^ 


0.134±0.022 


0.131±0.022 


0.108±0.018 


(2) 


PHENIX [2ii, 21] 


0.049±0.010 


0.047±0.010 


0.045±0.009 


V{2)l'^72) 


STAR [12, 2^ 


0.069±0.019 


0.067±0.019 


0.067±0.019 


(1)' (2) 


STAR [17. 181 


0.043±0.008 


0.043±0.008 


0.039±0.007 


A{2)/^(2) 


PHENIX [20, 21] 


0.072±0.017 


0.068±0.016 


0.074±0.017 


< > l\2) 


STAR [12, 24] 


0.039±0.011 






^l\2) 


STAR [12, 22] 


0.022±0.003 


0.021±0.004 


0.022±0.004 


^(2)/^(2) 


STAR [12, 28, 29] 


0.0093±0.0012 


0.0072±0.0011 


0.0060±0.0008 


^"A(2) 


STAR [12, 2^, 22] 


0.0014±0.0004 
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TT^^ : includes 50% of pions from decays of fi, S, S, A and their anti-particles 
K(2) • includes 100% of kaons from decay 
and 50% of kaons from Q and Cl decay 
: includes 100% of from cj) decay 
P(2) and p(2) : includes 100% of protons from the decay of A, S° and their anti-particles 
A(2) and A(2): includes 100% of the contribution from the decay of S, Q 
and their anti-particles 
ci-nd S^-j : includes 50% of the contribution from fl and fl decay 

As seen in Table I, in some cases the same ratios are at our disposal from various 
experiments. We include these multiply given ratios in our analysis as separate data 
points to be fit. 

V. RESULTS 

A. Centrality dependence of thermal parameters 

The global results of the fits are displayed in Fig. ^ As is evident from Fig. 
the chemical freeze-out temperature is roughly constant at around 165 MeV in Fit 
I. While the strangeness chemical potential is fairly flat at /i5 ~ 10 MeV (with 
10% variation, Fig. the baryon chemical potential (Fig. ^p) increases slightly 
with participant number in the range 33.4 MeV to 38.5 MeV (a 14% increase) in 
Fit I. Most striking, however, and the main result of our analysis, is the increase 
in strangeness saturation (Fig . H i) with increasing centrality for all fit types. This 
confirms previous arguing 12, |l5| . On a quantitative level, an analysis based on too 
few hadron species delivers non-reliable results, in particular also for 7^ (see Fig. di. 
Fit IV). Fits I and II, in contrast, deliver consistent results. Since we consider the 
phenomenologically introduced parameter 7^ as one possible indicator of deviations 
from equilibrium, the conclusion is that in central collisions equilibrium conditions 
for describing the chemical freeze-out are appropriate, while more peripheral colli- 
sions point to some off-equilibrium effects. 

To highlight the role of the parameter 7^ on the quality of the data fits, we exhibit 
in Fig. 121 the contours for fixed values of 7^ in the T — n plane. For definiteness, 



6 



_ 0.18r 
0) - 
£0.175 - 

I- 2 
0.17; 

0.165 ; 

0.16; 
0.155 
0.15; 

0.145 ; 
0.14 ; 

0.1 35 E 

_0.014r 

2.0.013 ; 

i = 

0.012 ; 

0.01 1; 
0.01 ; 

0.009 ; 
0.008 ; 
0.007 ; 
0.006 ; 

0.005 E 



100 150 200 250 300 




100 150 200 250 300 



^1.2p 
o 

U. 

c 

■i 1- 

g 

g:o.9- 

□ 

g 0.8- 
c 
o 
o> 

I 0.71 

3) 

0.6" 
0.5 c 



100 150 200 250 300 




100 150 200 250 300 



FIG. 1: The thermal parameters (a: temperature T, b: baryon chemical potential /jLb, c: 
strangeness chemical potential /i^, d: strangeness saturation factor 7^) as extracted from 
the experimental data for the various fit types (Fit I: solid, Fit II: dashed, Fit III: dotted. 
Fit IV: dash-dotted curves). 

only results for Fit I and peripheral collisions are shown. We have selected 7^ = 1 
(representing the assumption of full equilibrium), 0.85 (the optimum choice) and 0.5 
(as rather extreme choice), respectively. The corresponding x^/d.o.f are 24.1/12, 
17.2/11 and 62.0/12. The figure evidences that lowering 7^ noticeably shifts T up, 
while the up-shift of /i^ is less pronounced. The change in x^/d.o.f quantifies the 
relevance of the parameter 7^ for the fits. 



B. Detailed comparison with data 

After discussing the global trends, it is worthwhile to look into details. Instead of 
displaying the usual comparison of model results with data in one plot we consider 
each ratio individually to exhibit the centrality systematics. 
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FIG. 2: contours from Fit I of peripheral data. 7^ = 1 (0.5) in left (right) panel. In 
the middle panel, 7^ = 0.85 being the optimum choice. 

1. Particle- anti-particle ratios 

Fig. 121 summarizes such a comparison for the three centrality classes for the 
particle-anti-particle ratios. The data used from Table I are shown together with 
our model results. 

For a qualitative understanding one can use the Boltzmann approximation ig- 
noring widths. (At RHIC energies the use of Boltzmann statistics introduces errors 
in primordial yields at the level of 10% for pions, 1.3-1.5% for kaons and less than 
one percent for all other hadrons. Thus, given the intrinsic systematic uncertainties 
inherent in any thermal model analysis owing to the poorly constrained particle 
properties of the heavy resonances, quantum statistics is important only for pions.) 
The primordial density of hadron species i is then given by 

ni = j^^^s'^m^iTK2{mi/T) exp(/ii/T). (5.1) 

Thus, the ratio of particle to anti-particle primordial yields is simply, 

Tl ' 

- = exp(2/i,/T). (5.2) 

n-i 

Anti-particle-to-particle ratios therefore fix fiB,s/T. We note that these ratios are 
only slightly affected by feed-down (using the central Fit I best-fit parameters the 
percentage influence of feed-down is TT^^-j/vr^"^-): 1.0%, K^^^/ K^^y 2.9%, 0.5%, 
P{2)/Pi2y. 3.7%, A(i)/A(i): 0.7%, A(2)/A(2): 4.2%, E+^/El^y. 0.08%, and Cl/Q: 0%). 
Thus, the trends in these ratios translate directly into trends in /i^/T and fis/T, 
while they are roughly insensitive to 75 and T. 
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FIG. 3: Comparison between experiment and model results for particle-anti-particle ratios 

(a: ^(2)/^5)' ^- ^(t)/^(2)' P(i)Mi)' d: P(2)/P{2), e: A(i)/A(i), f: A(2)/A(2), g: ^'12)/ -{2)^ 
h: O^/r^^). Line codes as in Fig. ^ experimental results: STAR: stars, PHENIX: circles, 
PHOBOS: squares, BRAHMS: triangles, cf. Table I. 
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FIG. 4: As in Fig. 3 but for particle-to-pion ratios, (a: K,^)/^,^^ b: Kg/iT,^)^ c: p(i 



Let us now consider the individual ratios and begin with t^'{2)/'^'{2)- experi- 
mental ratio is consistent with being flat and slightly below or equal to unity, while 
our model yields a 7r^/7r+ ratio slightly above 1 (see Fig. Model and experi- 
ment, however, agree within errors. There is no significant difference between the 
results of the four fits. Setting /ig = drives the model ratio to very close to 1; 
deviations are due to feeding (since hb and ns are both greater than zero). With 
only strong decays included, the ratio is less than 1. Therefore it is the infiuence of 
the weak decays that pushes the ratio above 1. This ratio is insensitive to fis and 
fiB'- with a 5 MeV change in fiB the ratio varies by 0.2%, while a 5 MeV change 
in /i5 leads to a 0.3% change in the ratio. A 0.2 change in 75 furthermore affects 
the ratio by only 0.2%, while a 5 MeV change in T leads to only a 0.06% change in 
the ratio. Therefore, this ratio fixes nothing, and serves just as the motivation for 
setting /ig = 0. 

The ratio K'^.^/K^-j is above 1 (see Fig. Eb) and drives fis to positive values. 
Due to their much smaller errors, the STAR ratios are heavily weighted in the fit. 
This contributes to fis having a slight kink as seen in Fig. The effect of the 
above-quoted feeding accounts for the slightly different shape of the K^/K^ ratio 
compared with the fis/T trend which can be deduced from Fig. ^ With this ratio 
removed from the fit, fis/T would increase steeply with centrality, driven by the 
corrected A/A to double ratio. This gives H^/H^ an increasing trend. In 

order to keep the A/A ratio fiat, /is/T then similarly increases which causes p/p to 
drop. Therefore, the kaon ratio is what keeps fis fiattish and gives it the 'kink'. 

The p/p experimental ratios differ noticeably (see Fig. Efc and Fig. Ell), with 
STAR being the highest. These STAR data points are favored in fits since they have 
much reduced errors (consistently a factor of 3 less than those of PHENIX). These 
STAR ratios, however, are too high if the STAR A/ A ratio is to be simultaneously 
reproduced by the model (STAR p/p supports a lower /i^ as is seen when the A's 
are removed from the fit, as done in Fit IV). With A's excluded, the much smaller 
STAR errors lead to these p/p data points being strongly favored in the mid-central 
and peripheral fits. This leads to a strong increase of fis with centrality. Since the 
experimental A/A ratio (Fig. ^ and f) is fiattish, they support a similar behavior 
for Hb/T and Hs/T. 
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Excluding all but the STAR p/p points in the most central bin removes the strong 
increase in /i^ from the mid-central to central bins seen in Fit IV. Our suspicion that 
it is p/p which causes /i^ to increase is confirmed by repeating the fits with only anti- 
particle-to-particle ratios included (i.e., all particle-to-pion ratios excluded). Since 
these ratios are insensitive to 75, and T enters only together with the chemical 
potentials, in these reduced fits we fixed 75 = 1 and T = 165 MeV (both reasonable 
in light of the full results). The approximate trends in ^s/T and ^b/T observed in 
the full fits were reproduced using just these ratios. With the pjp ratio excluded, 
the steady increase in /i^ disappears; instead (ib/T decreases similarly to /xs/T as 
driven by the flatness in the A/ A ratio. 

The ratio (see Fig. 01;) depends strongly on fis/T and ^lb/T and is little 

affected by feed-down. Similar fractional errors to the uncorrected A/A ratio means 
both are treated equally in fits. 

Although consistent within errors, the model trend decreases slightly with central- 
ity, while the data are consistent with a weak increase. Comparing this experimental 
ratio with A/A, the S+/S^ to A/A double ratio increases. Thus, fis/T is driven up 
particularly from the mid-central to central bin. As mentioned earlier, the K^/K^ 
ratio leads to a kink in the trend for fis- This effect wins over the influence of the 
H to A double ratio which suggests an increase in fis/T. 

Unfortunately there is no centrality information on the ratio f2+/f2~. Such data 
would allow even better determination of fis- As seen in Fig. Eh, our model results 
agree fairly well with the data for the central bin. 



2. Particle-to-pion ratios 

The particle-to-pion ratios read in the Boltzmann approximation neglecting width 

n^ _ js,\ gi mf K2{mi/T) expjfij/T) 

i.e., qualitatively, mixed particle ratios allow the determination of 75 and T, and 
hence the /x's independently. Whereas the particle-to-anti-particle ratios feature 
just the combination fi/T, the primordial particle-to-pion ratios feature also T in 
the combination K2{mi/T)/ K2{mT^/T). This allows the temperature to be fixed 
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and then, given the values of fi/T fixed by the particle-to- anti-particle ratios, the 
chemical potentials are determined. However, feed-down plays an important role 
in these ratios, so that the final observed ratios differ from the ratios of primordial 
densities quite substantially. This makes a detailed analysis of the influence of these 
ratios on the parameters much more difficult than in the case of the particle-to- 
anti-particle ratios. In fact, for the best-fit parameters corresponding to Fit I of the 
central bin, more than 70% of tt+'s come from decays. 

A simple analysis shows how the ratio K2{mi/T)/ K2{m.,^/T) varies with temper- 
ature for a number of hadron species i: Lowering the temperature has the effect 
of bringing the particle-to-pion ratios down. Furthermore, the fractional effect of 
changing the temperature is greater for the heavier particles. 

After these general remarks let us consider in some detail the comparison of our 
model with data. The ratio K^^/tt^^) (Fig-Eti') is reasonably well reproduced in the 
fits. Since the errors on the PHENIX and STAR experimental ratios are similar, 
neither is heavily biased in the fit. The PHENIX data rise more rapidly than the 
STAR data (especially from the mid-central to central bin). Since the particle-to- 
anti-particle ratios support a centrality- independent fJ^s/T, the experimental /tt^ 
ratios drive an increase in strangeness saturation, with the PHENIX data supporting 
a steeper increase in 75 than the STAR data (see Fig. ^). In Fit I and II K~ /tt~ is 
over-predicted. In contrast, the uncorrected A/vr" and S^/tt^ are under-estimated 
(small fractional errors on S^/vr^ mean that this ratio strongly infiuences the fit). 
With the H's removed in Fit HI the STAR j-n^ data are well reproduced owing 
to the drop in T and 75. The big drop in this ratio when A's are excluded from 
the fit is due to the further drop in 75 and T in Fit IV seen in Fig. ^ The effect 
of this drop is somewhat lessened by the accompanying drop in /i^/T in the central 
and mid-central bins. Furthermore, the kaon to pion ratios are least affected by the 
temperature change in Fit IV. 

The increase in the ratio K^j-K^ with centrality (Fig. 0Jd) also drives 75 to in- 
crease. In the full fits (Fits I - HI) it is somewhat over-estimated (a similar result 
is found in |15|]), while the exclusion of the A's (Fit IV) results in much better 
agreement. This is owing to the reduced 75 in Fit IV. Comparing the results for 
the ratios j-n^ and K^/tt~ in Figs. and b one gets the impression that the 
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PHENIX K~ /t[~ data is compatible with STAR's K^/ti~ data, while the good fit 
in Fit 111 of star's K~ /t[~ data result in a ratio Kg/7i~ being at the upper limits 
of the error bars from STAR. The central < K*^ > /vr^ ratio reproduction by the 
model (Fig. 0^) follows the same trend as the K^. 

Let us now focus on the p/vr^ and A/vr^ ratios (Figs. Et, d, e and f). While the 
p/vr" ratio is over-estimated by the model in Fits 1 - 111, A/vr" is under-estimated 
except for the corrected A/ it" in the central bin in Fits 1 and 11. With A excluded, 
p/vr" is well reproduced, while A/vr" is greatly underpredicted. Therefore, p/vr" 
drives T down, while A/vr" tends to raise 7s and T. 

The flatness of the ratio (Fig- Eh) supports a centrality-independent 75. 

In Fits 1 and 11 there is fair agreement with the model, while this ratio is greatly 
under-estimated in Fit IV owing to the greatly reduced 75 and T. The increasing 
trend of the model results is due to the increase in 7s with centrality and the fact 
that ~ 7|. 

The increase of the ratio S^/tt^ (Fig. EJ) drives fis/T to increase as well as 75. 
On the other hand, it causes fis/T to drop. Again there is fair agreement with the 
model in Fits 1 and 11, while this ratio is greatly under-estimated in Fit IV. 

There is fair agreement of our model results with the fi^/vr" ratio (Fig. EI) in 
the central bin in Fits 1 and 11. However, with S and A removed agreement worsens. 

Our final conclusion is that Fits 1-11 give a quite satisfactory description of most 
ratios, with the exception of K^/tt^, < K*^ > /vr^ and corrected p/tt^. For the 
< K*^ > /tt^ and /n^ there are no data for mid-central and peripheral collisions. 
Thus our results may serve as predictions. 

C. Summary 

In summary, we present an analysis of the centrality dependence of thermal 
parameters deduced from hadron multiplicities within a statistical-thermal model. 
The grand-canonical formalism is employed with quantum statistics and resonance 
widths included for all of the fits. We rely entirely on published data for the reaction 
Au -|- Au at y/sNN = 130 GeV. Our aim is to employ the strangeness saturation 
factor 7s as one possible indicator for deviations from equilibrium conditions in de- 
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scribing the chemical freeze-out. Indeed, for non-central collisions the systematic 
deviation of 7^ from unity can be considered as a hint to off-equilibrium effects, 
while for central collisions and using a sufficiently complete set of available multi- 
plicity ratios at mid-rapidity there is no indication for deviations from equilibrium 
conditions. 

We present a detailed comparison of the model results with data. The trends 
in fiB,s/T are fixed by the anti-particle-to-particle ratios. The / K~ ratio is 
responsible for keeping /is/T essentially flat. This wins over the influence of the 
H to A double ratio which tends to increase fis/T. The p/p ratio drives ^b/T 
up, which wins over A/A's attempts to make fis/T and fis/T behave similarly. The 
temperature and strangeness saturation are fixed by the particle-to-pion ratios. The 
observed increase in 75 is driven predominantly by the increase in the /tt^ ratio 
and the flatness of fis/T as fixed by the K~^/K^ ratio. 

A possible extension of the present work is the analysis of saturation in the non- 
strange sector, parameterized by 7^ being a factor similar to 7^ 0]. Since our 
fits to the data deliver reasonably small values of we believe, however, that 
there is not too much room left for sizeable deviations of 7^ from unity. This 
issue deserves a separate investigation. More tempting is the analysis of data for 
y^sjvTV = 200 GeV. The explorative study in [l^ points to a similar behavior of 
7s as we find for ^/snn = 130 GeV. Even more interesting is to test heavy flavor 
production; due to their larger intrinsic masses, heavy flavors should have a higher 
sensitivity to the medium. In particular, hadrons with open and hidden charm may 
require a corresponding factor 7c. 
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